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REPORT  SUMMARY 

Broadband  data  are  accumulating  on  Berkeley  network  stations  for 
earthquakes  and  explosions.  The  variational  formulation  for  elastic 
surface  waves  has  been  extended  to  the  case  of  dissipative  media.  Broad¬ 
band  spectral  data  on  Sierra  Nevada  earthquakes  appear  to  be  unaffected 
by  Q,  yielding  data  crucial  to  discrimination  among  various  earthquake 
source  models.  The  Oroville  earthquakes  of  August  1975  have  been  in¬ 
vestigated.  The  main  shock  produced  ground  displacements  of  1-2  mm  at 
Jamestown  (JAS)  and  Whiskey town  (WDC)  broadband  stations,  providing  a 
test  of  system  performance  under  relatively  strong  motion.  Both  systems 
behaved  well,  producing  excellent  records  of  the  strong  motion  with  no 
indications  of  spurious  responses  commonly  seen  in  long  period  systems. 
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A  VARIATIONAL  FORMULATION  FOR  LOVE 
V.'AVES  IN  A  LAYERED  ANELASTIC  SOLID  * 


Walt  Silva 

Department  of  Geology  and  Geophysics 
University  of  California 
Berkeley,  California 


ABSTRACT 

A  method  is  presented  which  allows  the  variational  formulation 
for  elastic  surface  waves  to  be  extended  to  the  case  of  dissipative 
media.  With  this  formulation,  correct  to  second  order  in  the  loss, 
Rayleigh's  principle  can  be  applied  to  perturbations  of  the  Rayleigh 
quotient  to  yield  group  velocity  without  numerical  differentiation. 
Other  perturbations  can  be  used  to  find  the  change  in  phase  or  group 
velocity  due  to  changes  in  loss,  density,  or  moduli. 


*  In  press,  Geoph.  J.  R,  A.  S.,  1976, 


INTRODUCTION 


The  use  o£  Rayleigh's  principle  is  of  considerable  importance 
in  surface  wave  calculations.  As  suggested  by  Meissner  (1926)  and 
Jeffreys  (1959,  1961)  and  amended  by  Harkrider  (1968)  it  has  replaced 
numerical  methods  with  an  exact  formulation  for  calculating  group 
velocities.  In  addition,  Rayleigh's  principle  may  be  used  to  calcu¬ 
late  the  effect  on  the  velocity  dispersion  due  to  small  perturbations 
in  the  elastic  parameters.  This  type  of  information  is,  of  course, 
most  essential  in  solving  inverse  problems. 

However,  in  its  usual  form,  Rayleigh's  principle  is  inadequate 
for  dissipative  systems  since  its  use  requires  equating  the  time  average 
kinetic  and  potential  energies.  This  result  follows  from  the  virial 
theorem  and  imposes  a  vanishing  Lagrangian  for  the  system. 

The  purpose  of  the  following  development  is  to  construct  a 
Lagrangian  for  a  non-conservative  system  which  vanishes  and  then  apply 
Rayleigh’s  principle  to  the  resulting  equation.  This  can  be  done  by 
writing  a  Lagrangian  for  two  systems  operating  simultaneously  with  one 
losing  energy  as  the  other  gains  energy,  so  that  the  total  energy  is 
conserved.  An  estimate  of  the  error  introduced  in  combining  the  systems 
is  given  in  the  appendix. 

The  development  will  be  for  Love  waves  in  order  to  simplify  the 
equations.  The  extension  to  Rayleigh  waves  follows  naturally. 


FORMULATION 


As  a  working  example,  we  shall  consider  Love  waves  propagated 
along  the  surface  of  a  vertically  heterogenous  (layered)  attenuating 
half-space.  The  usual  equation  to  be  satisfied  in  each  layer  is 


A  V*  u  *■  yAr.  Va  U  -  yO  £C 

uu  / 


(1) 


where  u  is  the  transverse  particle  displacement,  p  is  the  medium 

are 

density,  and  p  and  p  /  the  conservative  and  non-conservative 
*  f  which 

Lame  parameters,  respectively,  /  are  in  general  frequency  dependent 
(Borcherdt  1971,  Silva  1975).  We  take  as  the  solution  to  equation  (1) 
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where  Px  is  the  spatial  frequency  such  that  the  horizontal  anelastic 
phase  velocity  c  is  given  by  c  =  ~  and  Ax  is  the  horizontal  at- 
tenuation  factor.  Application  of  the  usual  boundary  conditions  results 
in  the  continuity  of  both  Px  and  A^^  .  They  then  represent  eigenvalues 
to  be  determined  for  propagating  modes.  The  attenuation  factor  may  be 
defined  as 
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where  represents  the  effective  quality  factor  for  the  spatial 

decay  of  the  surface  wave.  Equation  (3)  can  then  be  Interpreted  as 
the  projection  of  the  attenuation  in  the  propagation  direction.  The 
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Direct  substitution  of  equation  (7)  into  equation  (8)  yields 
equations  (1)  and  (5).  We  now  substitute  the  solutions  (equations  (2) 
and  (5))  into  equation  (7)  and  then  space-time  integrate  the  Lagrangian 
density  to  yield  the  Lagrangian  of  the  combined  system.  In  the  time 
and  x  integrations,  where  periodic  solutions  were  assumed,  the 
integrations  result  simply  in  averages. 

Defining  L  as  the  Lagrangian  of  the  combined  system,  we  can 


write: 


L  -■  \ /° M  ) 


r 


/<n  V  \r+  J*. 

(9) 


~  J  rv*'  y. 

The  first  term  is  interpreted  as  a  time  average  kinetic  energy 
and  the  remaining  two  terms  as  a  time  average  potential  energy.  This 
allows  the  application  of  the  virial  theorem  (Moisewitsch,  1966) 
which  states  that  for  a  conservative  system  which  is  quadratic  in  its 
potential  energy  the  time  average  potential  and  kinetic  energies  are 
equal. 

We  may  then  put  equation  (9)  into  the  form  of  the  Rayleigh 


quotient: 


u'*  I-  =  (  >’/-  A‘)  I,  *  I 


with  th(_j  energy  integrals 


/>  V  V*  Jz  ;  xi  --  \  /•»  v  \r*  dz  ; 


XJ  \  /<*  Vr  '  1/“  *  *  clr 


Application  of  Rayleigh's  Principle 


Rayleigh's  principle  states  that  for  a  given  Rayleigh  quotient, 
as  in  equation  (10),  any  eigenvectors,  correct  to  first  order,  will 
yield  the  eigenvalue  appropriate  to  that  mode  correct  to  second  order. 
This  can  be  stated  more  precisely  by  considering  the  perturbation  to 
equation  (10)  due  to  a  snail  change  in  v  and  v+  . 

_  ,  fi4  <n> 

I 

where  the  perturbed  energy  integrals  are  given  by: 


f  (\r  v  *•)  ctz  ;  fr!  *  \  /*•  f(rr*)  ^ 


1  A  *  (V'  \r*‘)  Jz 


Making  use  of  equation  (11)  the  new  eigenvalues  correct  to  second  order 
may  be  calculated. 
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Using  equations  (10)  and  (11)  and  neglecting  second  order  in 
small  quantities  equation  (12)  becomes: 

U/  Uul0  r  (>,  f  Px  -  At9At)  lx  -  <13> 

and  using  c  =  ^  for  phase  velocity,  U  =  for  group  velocf^y, 

Fx  *  *x  ' 

we  can  rewrite  equation  (13)  in  the  following  form: 

U  --  .1  -  A%  9AX  JL*.  (iA) 

£  P*  J  c 

At  this  point  Rayleigh's  principle  may  again  be  invoked  and  the  elastic 
eigenvectors  may  be  used  in  the  calculation  of  the  energy  integrals. 

In  order  to  estimate  the  perturbation  in  AxffA^  due  to  a 

frequency  perturbation,  we  use  equation  (3)  and  neglect  the  term  contain- 

f  Q  —1 

ing  — .  This  is  correct  to  at  least  second  order  in  Qr  since  0T 
b  9  u  L  L 

is  a  slowly  varying  function  of  frequency  and  the  term  has  a  coefficient 


of  Q. 


(It  can,  in  fact,  be  estimated  using  a  difference  scheme 


since  Q  (w)  is  known  from  the  eigenvalues.) 

V/ith  this  approximation  equation  (14)  becomes 


l  +  Vi  *  9i 


lx 
C  I- 


which  reduces  to  the  usual  elastic  expression  as  ->  c£»  and  v  *►  v  . 
This  expression,  which  is  an  approximation  correct  to  at  least  second 
order  in  the  loss,  is  thought  to  be  preferable  to  numerical  differenti¬ 
ation  of  the  phase  dispersion  curve. 


In  order  to  calculate  the  effect  on  Q ^  due  to  a  small  change 
j  r.  Q  we  must  first  find  an  expression  for  )i  in  terms  of  y,.  and 


Qs  .  Thi.s  can  be  done  by  assuming  a  solution  of  equation  (1)  of  tin 


f  orn 


u  e 


-x  (  X-  x  -  ) 

V.  /V  A'  ' 


(16) 


where  -  i^A_  (Borcherdt,  1973).  Substituting  this  into 

equation  (1)  ,  and  defining  the  shear  velocity  vs  as  that  of  horaogeneou 
waves  (P  parallel  to  A)  we  can  write: 


/***  - 
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1  *  J 

By  substituting  equation  (17)  into  1^  and  I2  we  ate  able  to 
explicitly  calculate  the  change  in  the  Love  wave  quality  factor  S 
due  to  a  change  £QS  in  the  energy  integrals.  Vie  can  write  the 
entire  perturbed  Rayleigh  quotient  at  constant  frequency  as: 


(17) 
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where  we_have  defined 

v,3  \r  V*  f 


/>  Jx. 

2 


di 


J  +  CV 


-  j 


s<;  1.  *  e»g,v>»-f/y  »•*  1 

j  2  L  1  -  «v  *  J 

Using  equations  (3),  (10)  and  (11)  and  negle 


(19) 


dz. 


neglecting  second  order  terms 


in  small  quantities  we  can  reduce  equation  (18)  to: 
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(20) 


Writing  equation  (3)  as 

p*  f  fp,)  ,  f  r<?j 


we  can  express  the  perturbation  in  Ax  due  to  a  change  in 


iP.  f  (qA  *  P,  f  f  fft)  (22) 

Neglecting  the  first  term  since  it  is  at  least  second  order  in  the 
loss  we  then  have  from  equation  (20); 

p.*  -  a:  )  4  * 

(23) 

and  to  the  same  order  we  have  fron  equation  (3) 

£  Q,  £  -  2  Qta  i  f  (Ql  ) 

(24) 


£  £  (  Qi )  -* 

P* 


3  P, 


( 


This  formulation  then  allows  estimates  to  be  made  on  the  change 
in  the  Love  wave  quality  factor  or  attenuation  factor  due  to  a  small 
change  in  the  shear  qualify  factor.  Further  perturbations  involving 
vs  and  p  can  be  made  along  similar  lines  (Harkrider,  1968). 

The  extension  to  Rayleigh  waves  will,  of  course,  be  more  tedious  since 
equation  (9)  will  now  contain  four  integrals.  In  addition,  equation  (18)  now 
must  consider  Qp  as  well  as  Qs.  The  situation  is  unpleasant  though  tractable 
and  the  group  velocity  is  presently  being  programmed  by  the  author. 


11 


APPENDIX 


In  order  to  estimate  the  error  introduced  by  the  formulation 
in  this  paper, 

presented  /we  shall  consider  a  plane  shear  wave  propagating  along  the 
x  direction  in  a  Voigt  solid  (Kolsky,  1963). 

The  equation  of  motion  can  be  written 


h 


fw 

v  X 3 


A  '  '  II  r  />  ^ 
r  cl  X3  X 


with  the  solution. 


-  AX  ,  > 

W  =  e  Ocn  ( -  PX ) 


(i) 


(2) 


We  may  define  the  quality  factor  for  this  system  as: 
-  1 


0  "  r  ^r-  I  (3> 

I 

I 

Substituting  equations  (3)  and  (2)  into  (1)  we  arrive  at: 


P  -  A 


U/  /O 


/*  0  *  <?'a) 


(4) 


On  the  other  hand,  if  we  had  applied  the  Lagrangian  formulation  pre¬ 
in  this  paper 

sented/the  results  would  have  been: 


P  -  A 


Us*  <0 
/  — 


A 


(5) 


which  is  correct  to  second  order  in  the  loss. 
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P-  AND  SV-WAVE  CORNER  FREQUENCIES  OVER  LOW-LOSS  PATHS: 

A  DISCRIMINANT  FOR  EARTHQUAKE  SOURCE  THEORIES?* 

William  A.  Peppin 
Seismological  Laboratory 
University  of  Nevada 
Reno,  Nevada  89507 

and 

Gerald  W.  Simila 
Seismographic  Station 
University  of  California 
Berkeley,  California  94720 

ABSTRACT 

Broadband  data  recorded  at  Jamestown,  California,  yield  displacement 
spectrum  corner  frequencies  for  P-  and  SV-waves  from  18  trans-Sierra  Nevada 
events  in  the  range  3.2  <  M^  <  4.0.  The  ratio  R  of  P  to  SV  corner  frequency 
is  .96  ±  .17  Hz  for  these  events,  which  sample  a  wide  range  of  source  azi¬ 
muths  and  focal  mechanisms.  Assumptions  of  a  flat  far-field  source  spectrum 
and  frequency- independent  attenuation  in  the  1-20  Hz  band  lead  to  Qtf  not 
less  than  480.  Many  recent  source  theories  predict  factor-of-two  or  more 
variation  of  R  with  azimuth,  and  significant  deviation  of  its  average  value 
from  unity;  such  source  models  apparently  do  not  apply  to  earthquakes 

studied  here.  Previous  studies  of  R  by  Molnar  et  al.  (1973),  Stump  (1974), 
and  Bakun  et  al.  (1975)  reach  different  conclusions.  We  interpret  this 
discrepancy  not  as  a  contradiction,  but  as  an  indication  that  R  may  be  an 
extremely  useful  source-model  discriminant,  its  variation  indicating  that 
different  models  may  apply  in  different  source  regions.  Well-determined 
measurements  of  R  may  lead  to  a  clearer  understanding  of  the  physical 
processes  associated  with  earthquakes. 

*  This  is  a  shortened  version  of  an  article  with  the  same  title 
submitted  to  J.  Phs.  Earth,  1976. 
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INTRODUCTION 

Recent  papers  on  earthquake  source  theory  differ  quite  markedly  in 
their  predictions  of  the  frequency  content  of  P  and  S  waves  leaving  the 
source.  High-quality  data  have  been  used  by  Molnar  et  al.  (1973),  Stump 
(1974)  and  Bakun  et  al.  (1975)  to  investigate  the  frequency  content  of  P- 
and  S-waves,  but  with  differing  results.  The  present  study  addresses  these 
discrepancies  with  analysis  of  the  broadband  vertical-component  seismograph 
system  at  Jamestown  (JAS) ,  California.  The  system  is  described  in  previous 
Technical  Reports  for  this  Contract. 


DATA 

As  JAS  Is  situated  in  the  Sierra  Nevada  bathollth,  it  monitors  low- 
loss  granite  travel  paths  for  trans-Sierra  Nevada  earthquakes.  All  well- 
recorded  trans  Sierra  Nevada  events  from  06  July  1974  to  31  August,  1975, 
and  a  representative  sampling  of  the  many  Oroville  earthquake  aftershocks 
were  analyzed.  Epicenter  and  magnitude  information  is  given  in  Table  1 
Event  nos.  1-18).  The  nature  of  the  travel  paths  to  JAS  can  be  seen  in 
Figure  1,  which  also  shows  the  distribution  of  exposed  granites.  Repre¬ 
sentative  seismograms  subjected  to  analysis  are  shown  in  Figure  2. 


PROCESSING 

The  JAS  seismograms  were  digitized  from  FW  magnetic  tape  at  25  or  50 
samples  per  second  after  analog  alias  filtering.  A  10%  cosine  taper  was 
applied  and  spectral  moduli  were  computed  for  the  entire  vertical  Pg  or  Sg 
phase,  then  smoothed  (3-point  triangle  window).  The  instrument  response 
was  removed,  but  not  the  10-Hz  alias  filter  response.  Identical  processing 
was  applied  to  a  noise  sample  (different  for  P  than  for  SV)  preceding  each 


Figure  1^  Location  of  the  trans-Sierra  Nevada  events  (black  dots)  relative  to 
station  JAS  (black  star).  The  numbers  beside  each  dot  are  the  event  numbers 
in  Table  1.  The  shaded  areas  are  exposed  granitic  rocks  associated  with  the 
Sierra  Nevada  batholith.  Focal  mechanisms  representing  the  four  source  regions 
are  presented;  the  small  dots  within  each  solution  give  the  approximate  Pg 
travel  path  to  JAS. 
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1 


P  (left)  and  SV  (right)  seismograms  at  JAS.  The  bracket 
beneath  the  trace  is  2  sec  long. 


18 


event.  Example  spectra  are  shown  in  Figure  3.  Corner  frequencies  for 
all  spectra  were  estimated  by  overlaying  the  spectra  with  a  template,  flat 
at  long  period  and  decaying  to  high  frequency  as  frequency  squared  or  cubed 
(whichever  fit  the  individual  spectrum  better) .  The  black  dots  in  Figure  3 
show  the  corner  frequencies  selected.  These  are  given  together  with  an 
estimate  of  their  uncertainty  in  Table  1. 

RESULTS 

The  corner  frequencies  were  used  to  compute  R,  and  these  values  are  given 
in  Table  1.  The  observed  value  of  R  for  all  18  earthquakes  is  .96  ±  .17.  The 
earthquakes  represent  at  least  four  distinct  source  areas  at  widely  varying 
azimuths  from  JAS;  thus,  we  can  safely  say  that  these  values,  showing  equal 
P-  and  S-wave  corner  frequencies  to  a  precision  of  20%,  are  representative  of 
earthquakes  in  the  Sierra  Nevada  province.  The  values  of  R  are  plotted  in 
Figure  4  along  with  data  from  Molnar  et  al.  (1973)  and  Bakun  et  al.  (1975). 

The  main  point  of  this  paper  is  that  these  three  studies,  each  based  on  ex¬ 
cellent  data,  give  quite  different  and  distinctive  patterns  for  R  in  the 
three  source  regions  represented. 

DISCUSSION 

Propagation  Effects 

We  would  hope  the  observed  values  of  R  reflect  conditions  at  the  source, 
but,  since  the  travel  paths  are  long  (Table  1),  distortion  caused  by  propa¬ 
gation  might  be  severe.  We  consider  now  two  forms  of  this  distortion: 
reverberations  among  crustal  layers  and  anelastic  attenuation,  along  with 
evidence  suggesting  neither  effect  alters  observed  frequencies  appreciably, 
and  thus  our  observations  of  R  do  indeed  apply  to  the  source  characteristics 
of  these  earthquakes. 
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Figure  Summary  of  published  data  on  corner  frequencies  of  P-  and  S-waves. 
Small  round  dots  are  from  Molnar  jet  al.  (1973);  open  triangles  are  from  Bakun 
et  al .  (1975);  solid  black  squares  are  from  this  study.  Table  1. 


Layer  reverberation  would  appear  to  cause  little  distortion  of  the 
spectra  because  the  travel  paths  -  through  crystalline  rock  -  are  quite 
homogeneous.  This  is  illustrated  by  the  very  simple  form  of  the  P-phase 
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for  Events  3,  8,  9,  and  10  in  Figure  2. 

Anelastic  attenuation: 

For  frequency-independent  attenuation  in  the  1-20  Hz  band,  we  have 
spectral  correction  A(f)  (Asada  and  Takano,  1963) 

A(f)  «  exp{  -frft/Q}  (1) 

If  the  spectral  corner  frequency  fc  is  defined  as  that  frequency  at  which 
the  observed  spectrum  is  6  db  below  the  (assumed  flat)  long-period  spectral 
component,  we  can  estimate  Q  from  (1)  as 

Q  a  -Tifct/ln(.5)  =  4.5  fct  (2) 

provided  that  Q  is  responsible  for  the  corner  frequency.  Therefore,  if 
the  corner  frequencies  in  Table  1  were  caused  by  Q,  we  would  find  from  (2) 
that 

VQa  =  V6-  *  J*7  •  (3) 

Most  investigators  have  found  that  exceeds  by  a  factor  of  1.5  or 
more  (Takano,  1971;  Bakun  and  Bufe,  1975;  Kurita,  1975;  Bakun  et  al . ,  1975) 
This  discrepancy  plus  the  lack  of  correlation  of  R  with  epicentral  distance 
indicate  that  Q  is  not  causing  the  observed  corner  frequencies. 

[.over  bounds  for  Q  ,  based  on  (2),  are  presented  in  Table  1.  Values 

p 

are  quite  high  (500-1000)  compared  with  the  above-mentioned  determinations 

of  crustal  Q  in  western  California.  The  lowest  estimate  for  Q  in  Table  1 

p 

(Event  16)  is  probably  not  representative,  because  the  event  had  the  short¬ 
est  travel  path  (hence  lowest  value  of  t  in  (2))  of  the  earthquakes  studied 
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Close-In  observations: 

Peppln  (1975)  has  presented  close-in  (a  few  km)  spectra  of  Event  17, 
with  SV  and  SH  well-defined  corner  frequencies  of  10-12  Hz,  above  the  value 
given  in  Table  1.  This  discrepancy  may  be  a  source  effect  cause  by  rapid, 
local  sliding  on  the  fault  plane  (which  causes  the  high  close-in  corner  fre¬ 
quencies)  superimposed  upon  slower  sliding  over  a  larger  fault  segment  (which 
causes  the  lower  corner  frequency  seen  at  JAS) .  Such  a  mechanism  was  pro¬ 
posed  for  the  Oroville  main  shock  by  Hart  et  al.  (1975) . 

Coast-Range  observations: 

Spectral  corner  frequencies  apparently  unaffected  by  propagation  would 
be  unusual  for  such  long  crustal  travel  paths.  For  example,  Furuya's  (1969) 
observations  are  consistent  with  those  of  central  California  earthquakes  re¬ 
corded  ft  Berkeley  (more  high  frequency  evident  in  the  P-phase) .  We  have 
studied  the  Antioch  sequence  as  recorded  at  station  BRK,  35  km  SW,  a 
3-component  broadband  system  similar  to  that  at  JAS.  Three  events  provided 
excellent  spectral  data  with  higher  corner  frequencies  for  P  than  for  SV 
by  40%  (Table  1).  However,  Bakun  et  al.  (1975)  report  a  Q„  of  100  to  150 

p 

for  a  travel  path  through  the  Gabilan  granites,  which  represents  the  most 
attenuation-free  travel  path  known  in  west-central  California.  If  we 

insert  Q  =  150  in  (1) ,  we  see  that  the  3-Hz  S-wave  corner  frequences  for 

b 

these  events  are  very  probably  caused  by  attenuation.  Then,  if  >  Qg  » 
we  would  see  a  higher  corner  frequency  for  P  than  for  S  by  (1).  Indeed, 
because  of  (1)  and  because  generally  Q  >  Q  ,  we  must  beware  of  interpreting 
corner  frequencies  in  terms  of  source  parameters  when  the  P-wave  gives  a 
higher  corner  frequency  than  S;  such  an  effect  almost  always  can  be  ex¬ 


plained  by  attenuation. 
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Implications  for  Earthquake  Source  Theory: 

We  assume  propagation  effects  have  not  significantly  altered  the 
observed  corner  frequencies  presented  in  Table  1,  so  that  these  corner 
frequencies  measure  directly  some  source  parameters.  In  all  recently- 
published  source  theories,  the  dimension  of  faulting  controls  the  P-  and 
S-wave  corner  frequencies,  and  in  all  but  two  the  predicted  variation  with 
azimuth  amounts  to  a  factor  of  2  or  more.  Thus,  the  resolution  of  the  data 
in  Table  1  is  adequate  to  discriminate  between  source  models  (subject  to 
the  assumption  above).  Berckhemer  and  Jacob  (1968)  employ  subsonic  rupture 
propagation  velocity,  but  concentration  of  energy  from  the  center  of  a  cir¬ 
cular  dislocation  to  find  equal  P-  and  S-wave  corner  frequencies.  Brune 
(1970,  1971)  and  Hanks  and  Wyss  (1972)  employ  an  infinite  rupture  speed  on 
a  circular  fault  to  find  higher  corner  frequencies  for  P  than  for  S  by  the 
ratio  of  the  P-  and  S-wave  velocities.  Molnar  et  al.  (1973)  use  uniform 
rupture  at  high  propagation  velocity  on  a  circular  dislocation  to  find 
higher  values  for  P-  than  for  S-wave  corner  frequencies.  Sato  and  Hirasawa 
(1973)  use  subsonic  rupture  propagation  on  a  circular  dislocation  with 
greater  slip  near  the  center  to  find  higher  (>  50%)  corner  frequencies  for 
P  than  for  S.  Burridge  (1975)  considers  self-similar  rupture  on  a  circular 
fault  lacking  cohesion  at  the  P-wave  velocity.  His  formula  (59)  for  the 
ratio  of  the  P-  to  SV-phase  corner  frequency  is  1.2  averaged  over  the  focal 
sphere,  but  the  ratio  is  2  or  more  on  50%  of  the  focal  sphere.  Madariaga 
(1975)  solves  the  problem  of  slip  on  a  circular  crack  when  stress  drop  is 
specified.  The  slip  is  concentrated  near  the  center  of  the  crack,  and 
gives  slightly  higher  (10-20%)  P-corner  frequencies  for  subsonic  rupture 
propagation.  In  contrast.  Savage  (1972)  uses  a  Haskell 
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fault  model  with  uniform,  bilateral,  subsonic  rupture  propagation  at  90%  of 
the  S-wave  speed  and  uniform  fault  slip;  he  finds  higher  primary  corner  fre¬ 
quencies  for  S  than  for  P  by  a  ratio  (averaged  over  the  focal  sphere)  of  1.65 
to  1.  Similarly,  Dahlen  (1974)  uses  a  combined  dislocation  and  self-similar 
crack  model  at  subsonic  rupture  propagation  with  uniform  fault  slip;  he  finds 
higher  corner  frequencies  for  S  than  for  P  (averaged  over  the  focal  sphere) 
by  the  ratio  2.38  to  1  or  1.75  to  1,  depending  on  whether  the  rupture  propa¬ 
gation  velocity  is  90%  or  80%  of  the  S-Wd,^  velocity.  In  summary,  models  which 
have  uniform  slip  and  subsonic  rupture  propagation  on  circular  or  bilateral 
faults  predict  higher  corner  frequencies  for  S  than  for  P  at  all  azimuths, 
while  models  which  employ  supersonic  rupture  propagation  velocity  or  slip  con¬ 
centrated  near  the  fault  center  predict  higher  corner  frequencies  for  P  than 
for  S  at  all  azimuths  (on  the  average);  Table  2.  Of  the  above-mentioned  models, 
only  those  of  Berckhemer  and  Jacob  (1968)  and  Madariaga  (1975)  cannot  be  ex¬ 
cluded  by  the  data  of  Table  1;  but  note  that  the  former  authors  employ  what 
may  be  an  unrealistic  assumption,  namely  rupture  propagation  velocity  decreasing 
as  the  faulting  front  moves  out. 

The  quantity  R,  then,  appears  to  be  a  rather  sensitive  source-model  dis¬ 
criminant.  Our  data  could  be  explained  by  the  models  of  Berckhemer  and  Jacob 
(1968)  or  Madariaga  (1975),  but  not,  for  example,  by  the  model  of  Savage  (1972). 
Bakun  ejt  al.  (1975)  find  the  opposite:  Savage's  model  can  describe  their  data 
quite  well.  None  of  these  three  seems  capable  of  explaining  the  high  values 
of  R  found  by  Molnar  et^  al.  (1975),  but  Brune's  (1970,  1971)  model  with  a  modi¬ 
fication  by  Hanks  and  Wyss  (1972)  can.  We  saw  in  Table  2  that  R  depends 
largely  upon  the  rupture  propagation  velocity  on  the  fault;  it  is  physically 
reasonable  that  this  quantity  might  vary  in  different  source  regions  (e.g. 
slower  in  areas  of  great  heterogeneity,  faster  on  smooth,  well-developed  fault 
traces  like  the  San  Andreas).  That  is,  we  should  not  be  surprised  to  find  that 
no  one  model  can  accurately  predict  the  totality  of  earthquake  phenomena  (but 


note  in  Table  2  the  features  common  to  the  models:  each  can  be  formulated  as 
a  specific  case  of  Haskell  dislocation  theory). 

The  results  of  this  investigation  can  be  explained  by  another  model.  If 
the  corner  frequency  is  determined  by  the  time  history  of  motion  at  the  source 
for  these  events,  then  equal  and  non-varying  corner  frequencies  for  P  and  S 
would  result  at  all  azimuths.  This  particular  model  seems  physically  reasonable 
for  small-to-moderate  earthquakes;  it  implies  a  source  time  duration  that  is 
long  compared  with  the  rupture  propagation  time  acorss  the  source.  This  ex¬ 
planation  was  considered  by  Johnson  and  McEvilly  (1974)  for  Central  California 
earthquakes,  and  would  be  consistent  with  the  findings  of  Archuleta  and  Brune 
(1975)  based  on  model  experiments  in  foam  rubber.  More  complete  azimuthal 
coverage  would  be  needed  to  test  such  a  model. 

CONCLUSIONS 

P-  and  SV-wave  spectral  corner  frequencies  for  18  trans-Sierra  earthquakes 
give  a  ratio  R  of  P-wave  to  S-wave  corner  frequency  of  .96  -  .17.  These  values 
appear  to  be  independent  of  propagation  effects,  and  thus  are  indicative  of 
conditions  at  the  source  of  the  earthquakes.  The  SV-wave  spectra  permit  a 
minimum  estimate  of  Qg  of  480,  which  is  high  compared  with  other  California 
crustal  paths  previously  studied.  The  precision  with  which  we  have  determined 
R  allows  us  to  eliminate  all  but  two  of  the  recently  published  source  models 
as  applicable  to  these  earthquakes;  in  this  way  we  find  that  R  is  a  successful 
discriminant  for  competing  source  models. 
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TABLE  2.  Essentials  of  Recent  Source  Theories 
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ABSTRACT 

The  characteristics  of  the  Oroville,  California  earthquake  sequence 
of  1975  are  presented.  Historically,  the  Oroville  area  is  one  of  low 
seismicity;  the  largest  earthquake  in  the  region  occurred  with  a  magnitude 
5.7  in  1940  some  50  km  north  of  Oroville.  This  low  seismicity  continued 
through  initial  loading  of  nearby  Oroville  Dam  in  1967-68  and  until  the 
first  foreshock  on  June  28,  1975.  Twenty-one  foreshocks  (M^  _>  1.6),  the 
largest  of  magnitude  4.7,  preceded  the  magnitude  5.7  main  shock  of  01  August. 
All  foreshocks  and  aftershocks  of  3.0  were  located  using  seismographs 

operated  by  the  University  of  California  at  Berkeley,  USGS,  and  the 
California  Department  of  Water  Resources.  The  aftershock  region  covers 
an  area  approximately  14  km  by  10  km  southeast  of  the  city  of  Oroville. 

The  depth  distribution  of  the  earthquakes  indicates  a  west  dipping  fault 
plane.  The  b  value  of  0.61  shows  the  sequence  to  be  rich  in  larger  mag¬ 
nitude  aftershocks.  Similar  b  values  have  been  determined  for  other 
aftershock  sequences  in  California,  such  as  a  sequence  near  Coalinga  in 
August  1975.  The  aftershock  occurrence  rate  follows  an  Omori  relation 
with  n(t)<*t  Apparent  variability  in  the  earthquake  mechanisms  of 

the  series  makes  interpretation  of  composite  fault  plane  solutions  difficult, 
but  the  data  indicate  normal  faulting  striking  NNW  and  a  west  dip  of  about 
30°  -  40°,  the  Sierra  Nevada  moving  up  with  respect  to  the  Great  Valley. 

*  In  press,  Bull.  Seism.  Soc  Amer.,  Aug.,  1976. 
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On  August  1,  1975  at  2020  GMT(1320  PDT)  a  magnitude  5.7  earthquake 
occurred  12  kilometers  southwest  of  Oroville  Dam  (Figure  1)  and  caused 
some  damage  in  the  city  of  Oroville  and  vicinity.  This  earthquake  was 
preceded  by  foreshocks  beginning  on  June  28,  1975,  the  largest  of  which 
had  Ml  =  3.5  (Table  1);  the  sequence  appeared  to  have  died  out  (only 
five  events  between  July  8  and  August  1),  when  on  August  1,  11  foreshocks 
with  magnitude  less  than  4.8  occurred  within  five  hours  of  the  main  shock. 

The  results  of  seismic  monitoring  near  Oroville  are  shown  in  Figure  2. 
The  histogram  shows  numbers  of  earthquakes  within  10  kilometer  increments 
of  radius  out  to  50  kilometers  from  ORV.  The  increase  in  seismicity  in 
Figure  2  from  1968  to  1975  reflects  the  5-6  fold  increase  in  magnification 
achieved  with  the  switch  from  paper  (U.  C.  Berkeley)  to  16  mm  film  (DWR) 
records. 

The  elevation  of  the  water  surface  of  Lake  Oroville  is  also  shown  in 
Figure  2.  From  November  1967  until  March  1968,  the  reservoir  filled  to 
half  its  capacity  (water  surface  elevation  229  m) .  The  reservoir  was  first 
filled  to  capacity  in  July  1969.  The  largest  annual  level  fluctuations 
occurred  in  the  winter  of  1974-75  when  the  reservoir  level  reached  a  low 
of  230  m  elevation.  This  draw-down  is  equivalent  to  48%  of  the  reservoir 
capacity.  In  prior  years  the  draw-down  was  typically  about  25%.  During 
the  seven  year  period  of  large  water  storage,  the  histogram  indicates  no 
significant  fluctuations  in  the  local  seismicity  near  to  the  dam  (<  30  km). 
The  very  low  level  of  local  activity  continued  through  the  1974-75  rapid 
draw  down  and  recharge. 


On  June  28,  the  foreshocks  of  the  1975  Oroville  series  began.  These 
events  are  shown  on  the  histogram  of  Figure  2  in  the  0-10  and  10-20  km 
ranges.  As  a  reaction  to  the  earthquake  of  June  28,  the  portable  seis¬ 
mograph  stations  Butte  City  (BUT)  and  Sutter  Butte  (SUT)  were  established 
on  July  1,  1975.  Locations  are  shown  in  Figure  3. 

The  magnitude  of  the  mainshock  is  estimated  to  be  5.7  on  the  Richter 
scale.  This  value  was  obtained  by  averaging  the  individual  estimates  from 
Wood-Anderson  seismographs  at  Areata,  5.7;  Berkeley,  5.7;  Stanford,  5.7 
(T.  Hanks  personal  communication);  and  Mount  Hamilton,  5.6  (one  component). 
No  estimate  was  obtained  from  the  instrument  at  Mineral  as  its  recordings 
were  off-scale. 

The  hypocentral  parameters  and  Richter  magnitudes  of  the  51  earth¬ 
quakes  (Ml  ^3.0)  of  the  Oroville  sequence  occurring  between  June  28  and 
August  16,  1975  are  listed  in  Table  1.  Locations  of  these  earthquakes 
are  shown  in  Figure  4.  The  bars  in  the  figure  represent  standard  errors. 

The  aftershock  region  covers  an  area  approximately  10  km  by  14  km  centered 
8  km  southeast  of  Oroville.  The  earthquakes  are  generally  deeper  in  the 
western  part  of  the  aftershock  zone,  with  depths  ranging  from  1.5  to  9.2  km. 

The  Oroville  sequence  through  August  11  consisted  of  21  foreshocks, 
the  main  shock,  and  292  aftershocks  with  Richter  magnitude  of  1.6  or  above. 
Twelve  aftershocks  with  Ml  >  4.0  were  recorded  during  the  same  time  period. 

A  graphical  representation  of  the  earthquake  sequence  is  shown  in  Figure  5. 

The  frequency  of  occurrence  shown  in  Figure  6  may  be  fitted  to  Omori's 
relation:  n(t)  *  A/(t+C)^  .  Aftershocks  were  divided  into  cells,  each 


containing  15  earthquakes,  and  the  average  occurrence  time  plotted  for  each 
cell.  The  least  squares  solution  gave  p  -  0.70. 


Magnitude  occurrence  data  are  summarized  in  Figure  7  where  N  represents 


the  cumulative  number  of  earthquakes  of  magnitude  greater  than  or  equal  to 
M^.  The  Ml  values  are  plotted  in  magnitude  increments  of  one  tenth.  The 
equation  log  N  s  a-bMk  was  fitted  by  least-squares  and  yielded  b  ■  0.61. 

The  main  event  magnitude  is  included  in  the  solution.  The  b  value  changes 
by  less  than  2%  if  the  main  event  is  removed  or  if  only  events  of  Ml  >  2.8 
are  used. 

Composite  fault  plane  solutions  were  determined  for  earthquakes  in  the 
northern,  west  central,  and  southern  parts  of  the  aftershock  zone.  Due  to 
the  limited  data,  well-resolved  fault  planes  were  difficult  to  obtain.  The 
northern  and  west  central  parts  of  the  aftershock  zone  clearly  indicate 
normal  faulting  with  the  mountains  to  the  east  up  with  respect  to  the  valley 
to  the  west.  The  group  of  4  earthquakes  to  the  north  (Figure  8a)  yields  a 
solution  with  a  strike  of  between  N  10°  W  to  N  20°  E  and  a  west  dip  of 
approximately  40°.  The  west  central  section  (Figure  8b)  indicates  a  fault 
strike  range  of  N  10°  W  to  N  30°  W  with  a  west  dip  of  30°.  No  coherent 
solution  could  be  found  for  the  southern  group. 

Testing  JAS  and  WDC  broadband  systems  for  large  motions: 

The  Oroville  shock  provided  an  excellent  test  of  the  broadband  flat 
displacement  channels  at  JAS  (195.5  km,  151.3°  Az)  and  WDC  (149.6  km, 

324.4°  Az) .  Figures  8  and  9  show  magnetic  tape  playouts  from  the  stations. 
High-  and  low-gain  displacement  channels  differ  in  gain  by  200,  as  indicated. 

The  low-gain  channels  behaved  very  well,  recording  nearly  2  mm  of 
ground  displacement  with  no  indication  of  spurious  seismometer  response 
such  as  spring  resonance  or  impulsive  responses  to  zero  position  shifts 
frequently  seen  in  long  period  systems. 
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The  high-gain  displacement  channels  also  behaved  well.  Interestingly, 
the  long  period  (10-20  sec)  motion  generated  by  this  shock  was  visible  on 
these  records  for  several  hours,  possibly  reflecting  surface  wave  energy 
trapped  and  resonating  within  the  Great  Valley. 

In  conjunction  with  excellent  data  obtained  at  short  distances  by 
local  networks,  the  regional  broadband  recordings  will  provide  excellent 
information  on  this  earthquake. 
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Earthquake  data  for  events  with  ML13.0 


Date 

Origin 

Time 

(OT) 

0(T)sec 

Magni¬ 

tude 

ml 

Latitude 

_ 

a(cf>)km 

Longitude 

X 

o(X)km 

Depth 

h 

0(h)  km 

1 

28 

JUN 

75 

041953.72 

0.13 

3.5 

39< 

'28.33’ 

1.19 

121' 

‘31.45’ 

1.05 

7.56km 

1.59 

2 

01 

AUG 

75 

154537.81 

0.07 

3.8 

39 

26.98 

0.67 

121 

31.87 

0.67 

7.27 

1.06 

3 

01 

AUG 

75 

162717.81 

0.14 

4.7 

39 

26.29 

1.19 

121 

32.25 

1.20 

4.89 

2.17 

4 

01 

AUG 

75 

172650.12 

0.09 

3.0 

39 

27.71 

0.85 

121 

32.27 

0.83 

8.57 

1.22 

5 

01 

AUG 

75 

202004.75 

0.08 

4.5 

39 

26.33 

0.76 

121 

31.71 

0.76 

8.01 

1.17 

6 

01 

AUG 

75 

202012.85 

5.7 

39 

26.33 

121 

31.71 

7 

01 

AUG 

75 

2025 

4.7 

39 

26.33 

121 

31.71 

8 

01 

AUG 

75 

2029 

4.6 

39 

26.33 

121 

31.71 

9 

01 

AUG 

75 

203239.84 

0.15 

3.0 

39 

26.71 

1.62 

121 

30.40 

1.44 

4.82 

3.21 

10 

01 

AUG 

75 

2037 

3.5 

39 

26.71 

121 

30.40 

11 

01 

AUG 

75 

2045 

3.0 

39 

26.71 

121 

30.40 

12 

01 

AUG 

75 

204618.40 

0.14 

3.8 

39 

28.40 

1.56 

121 

30.02 

1.59 

6.13 

2.69 

13 

01 

AUG 

75 

210539.84 

0.04 

3.0 

39 

25.98 

0.38 

121 

29.25 

0.43 

6.74 

0.63 

14 

01 

AUG 

75 

212150.65 

0.07 

4.1 

39 

26.53 

0.64 

121 

31.70 

0.64 

7.76 

0.99 

15 

01 

AUG 

75 

212559.02 

0.04 

3.3 

39 

28.44 

0.42 

121 

31.08 

0.41 

7.08 

0.63 

16 

01 

AUG 

75 

212924.12 

0.08 

3.6 

39 

27.12 

0.69 

121 

32.92 

0.80 

6.55 

1.21 

17 

01 

AUG 

75 

234440.98 

0.03 

3.4 

39 

29.18 

0.25 

121 

31.34 

0.24 

7.51 

0.35 

18 

02 

AUG 

75 

005248.47 

0.03 

3.8 

39 

29.06 

0.25 

121 

30.55 

0.24 

7.27 

0.36 

19 

02 

AUG 

75 

063157.19 

0.06 

3.2 

39 

26.84 

0.50 

121 

29.01 

0.38 

5.68 

0.79 

20 

02 

AUG 

75 

101153.68 

0.04 

3.1 

39 

29.40 

0.37 

121 

30.70 

0.35 

7.22 

0.52 

21 

02 

AUG 

75 

104900.12 

0.07 

3.3 

39 

25.66 

0.52 

121 

28.41 

0.43 

5.54 

0.85 

22 

02 

AUG 

75 

115150.74 

0.12 

3.4 

39 

28.36 

0.89 

121 

29.30 

0.91 

2.49 

1.86 

23 

02 

AUG 

75 

144438.71 

0.04 

3.2 

39 

25.03 

0.32 

121 

29.38 

0.33 

5.00 

0.59 

24 

02 

AUG 

75 

165145.11 

0.05 

3.7 

39 

25.07 

0.43 

121 

29.09 

0.44 

6.48 

0.75 

25 

02 

AUG 

75 

172429.23 

0.07 

4.3 

39 

28.42 

0.60 

121 

28.25 

0.59 

5.84 

0.98 

26 

02 

AUG 

75 

174324.13 

0.06 

4.0 

39 

28.66 

0.51 

121 

28.41 

0.49 

5.52 

0.83 

27 

02 

AUG 

75 

195836.85 

0.05 

3.1 

39 

26.89 

0.44 

121 

32.20 

0.44 

7.10 

0.70 

28 

02 

AUG 

75 

202216.32 

0.05 

5.1 

39 

26.69 

0.35 

121 

27.76 

0.36 

3.62 

0.71 

29 

02 

AUG 

75 

203548.55 

0.06 

3.9 

39 

28.26 

0.52 

121 

28.89 

0.51 

5.98 

0.84 

30 

02 

AUG 

75 

205855.70 

0.06 

3.8 

39 

25.93 

0.50 

121 

27.97 

0.51 

5.71 

0.90 

31 

02 

AUG 

75 

2059 

5.2 

39 

25.93 

121 

27.97 

32 

02 

AUG 

75 

214001.34 

0.07 

3.9 

39 

25.62 

0.57 

121 

28.39 

0.58 

5.39 

1.04 

1.04 
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03 

AUG 

75 

010305.77 

0.04 

4.6 

39 

29.26 

0.40 

121 

31.06 

0.38 

7.98 

0.55 

34 

03 

AUG 

75 

024708.80 

0.04 

4.1 

39 

28.68 

0.33 

121 

30.06 

0.32 

6.79 

0.50 

35 

04 

AUG 

75 

094745.01 

0.06 

3.5 

39 

25.27 

0.61 

121 

31.36 

0.45 

7.57 

0.88 

36 

05 

AUG 

75 

022857.35 

0.05 

3.3 

39 

24.91 

0.39 

121 

29.46 

0.33 

7.17 

0.63 

37 

05 

AUG 

75 

204424.50 

0.17 

3.2 

39 

24.87 

1.35 

121 

30.87 

1.15 

6.79 

2.23 

38 

06 

AUG 

75 

035029.94 

0.20 

4.7 

39 

28.73 

1.94 

121 

31.46 

1.73 

7.69 

2.51 

39 

06 

AUG 

75 

162547.93 

0.14 

3.1 

39 

26.73 

1.24 

121 

27.58 

1.14 

8.49 

1.30 

40 

06 

AUG 

75 

164152.08 

0.02 

3.6 

39 

29.80 

0.17 

121 

31.76 

0.17 

8.40 

0.16 

41 

06 

AUG 

75 

210033.46 

0.11 

3.0 

39 

26.37 

0.95 

121 

29.10 

0.89 

9.20 

0.98 

42 

06 

AUG 

75 

2100 

3.0 

39 

26.37 

121 

29.10 

43 

07 

AUG 

75 

203120.38 

0.15 

3.1 

39 

30.99 

1.56 

121 

31.98 

1.48 

8.68 

1.71 

44 

08 

AUG 

75 

070050.09 

0.05 

4.9 

39 

30.13 

0.51 

121 

30.72 

0.46 

7.64 

0.62 

45 

08 

AUG 

75 

133753.87 

0.08 

3.2 

39 

29.79 

0.70 

121 

29.37 

0.64 

6.29 

0.68 

46 

08 

AUG 

75 

190327.23 

0.09 

3.1 

39 

23.60 

0.65 

121 

29.60 

0.62 

5.60 

0.91 

47 

09 

AUG 

75 

073847.48 

0.10 

3.0 

39 

24.80 

0.74 

121 

29.00 

0.70 

7.46 

0.90 

48 

11 

AUG 

75 

024016.73 

0.35 

3.0 

39 

27.72 

2.28 

121 

26.53 

1.85 

1.52 

3.89 

49 

11 

AUG 

75 

061136.34 

0.12 

4.3 

39 

26.78 

0.91 

121 

28.85 

0.81 

4.26 

1.24 

50 

11 

AUG 

75 

155905.26 

0.19 

3.6 

39 

28.22 

1.40 

121 

33.23 

1.39 

6.39 

1.68 

51 

16 

AUG 

75 

054809.38 

0.03 

4.0 

39 
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FIGURK  1:  Station  locations  and  historical  seismicity  within  60  km 


of  Oroville. 


FIGURE  2:  Local  earthquakes  in  the  Oroville  area  and  Lake  Oroville 
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FIGURE  4:  Earthquake  locations  of  the  Oroville  sequence.  Standard 


errors  in  location  estimated  by  least  squares  are  indicated 
by  bars  through  the  epicenters.  Dashed  diamonds  represent 
earthquakes  whose  locations  are  assumed  to  be  those  of  their 
immediate  foreshocks.  The  dashed  lines  indicate  the  lines 


of  surface  cracks  observed  after  the  main  shock 
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6:  Logarithmic  representation  o£  Figure  5  fitting  an  Omori  type 
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FIGURE  7:  Plot  of  log  N  vs  ML  (ML  >  1.6).  The  solid  line  is  the  least 


squares  fit  to  the  data. 


Figure  8.  Broadband  recordings  of  the  01Aug75  Oroviile  earthquake  at  JAS,  vertical 
component.  Low-gain  (upper)  channel  is  +2mm  full  scale  ground  motion; 
high-gain  channel  is  +J0microns  full  scale.  Arrows  indicate  P-wave 
arrivals  for  foreshock  and  main  shock. 
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